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ABSTRACT 

We describe a new search for diffuse ultrahigh energy gamma-ray emission associated with 
molecular clouds in the galactic disk. The Chicago Air Shower Array (CAS A), operating in 
coincidence with the Michigan muon array (MIA), has recorded over 2.2 x 10 9 air showers from 
April 4, 1990 to October 7, 1995. We search for gamma rays based upon the muon content of 
air showers arriving from the direction of the galactic plane. We find no significant evidence for 
diffuse gamma-ray emission, and we set an upper limit on the ratio of gamma rays to normal 
hadronic cosmic rays at less than 2.4xl0~ 5 at 310 TeV (90% confidence limit) from the galactic 
plane region: (50° < l\\ < 200°; —5° < b\\ < 5°). This limit places a strong constraint on models 
for emission from molecular clouds in the galaxy. We rule out significant spectral hardening in 
the outer galaxy, and conclude that emission from the plane at these energies is likely to be 
dominated by the decay of neutral pions resulting from cosmic rays interactions with passive 
target gas molecules. 



Subject headings: gamma ray: observations — cosmic rays — ISM: clouds — Galaxy: 
fundamental parameters 



Introduction 



Diffuse gamma rays arriving from regions of enhanced density, such as giant molecular clouds, are 
powerful tracers of the interactions of cosmic rays with matter in the galaxy. We expect diffuse emission 
from molecular cloud regions where the permeating flux of high energy cosmic rays interacts with the 
nuclei of gas atoms (mostly hydrogen) to produce new hadronic particles, including neutral pions, which 
subsequently decay into gamma rays. In this manner, the gas in molecular clouds acts as a passive target, 
converting some fraction of impinging cosmic rays into gamma rays. Such diffuse emission has been recorded 
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by several space-borne detectors including the EGRET experiment aboard NASA's Compton Gamma Ray 
Observatory - indeed the all-sky flux at 100 MeV is largely (~ 90%) diffuse emission from molecular clouds 
in the galactic plane. Most of this emission is contained within a narrow band along the galactic equator 
(fen < ±5°), and is sharply concentrated towards the galactic center. The EGRET results have been used, 
together with radio data, to develop a three-dimensional model of both the gas and cosmic ray densities 



in the galaxy (Bcrtsch et al. 1993, Hunter et al. 1997). This model, which estimates the contributions 
to diffuse emission from electron bremsstrahlung, nucleon-nucleon interactions, and inverse Compton 
processes, accurately matches the observed emission in detail as seen by EGRET for all galactic longitudes 
over the energy range from 30 MeV to about 1 GeV. (The model deviates somewhat from observations at 



energies above 1 GeV, possibly due to a change in the cosmic ray spectrum; see Mori 1997 ). The EGRET 
results support a galactic origin of cosmic rays, presumably as a result of shock acceleration in supernova 
remnants (SNR). 

Much less is known about the nature and distribution of cosmic rays at energies above 100 TeV 
(10 14 eV). Past efforts have mostly focused on searches for point sources of gamma rays using large-area 
ground-based air shower detectors. However, after many sensitive searches, there exists no compelling 
evidence for the detection of even a single persistent gamma-ray point source at these energies (see e.g., 
Alexandreas et al. 1991, Cronin et al. 1993, Borionc ct al. 1997a, Borionc et al. 1997b| ). 



We can estimate the ultrahigh energy diffuse flux expected from clouds in the galaxy by approximating 
the cosmic ray density and the matter density within any particular molecular cloud as uniform. In this 
case, the diffuse gamma-ray flux from this interaction can be expressed as: 



■M> E) = i-g 7 (> E)> < N H >, 



(1) 



where < Nh > is the total hydrogen column density in a given direction through the cloud and q 7 (> E) is 
the source function for gamma-ray emissivity: 



q 7 {> E) = Am J OR (> E) ■ T{a { 



(2) 



Here J C r(> E) is the ambient integral flux of normal hadronic cosmic rays and T is a function of the total 
inelastic cross section, (Tj ne j, the integral spectral index of cosmic rays, 7, and a correction factor, j 'a which 
takes into account that some primaries and targets are nuclei, and not protons. The function T varies 
slowly with energy, and is largely independent of the form of the cosmic ray spectrum at ultrahigh energies, 
(see e.g., Dcrmcr 1986, |Gaisser and Stancv 1991 , Aharonian 1991, Bcrczinsky and Kudryavtsev 1990 for 
details concerning diffuse emissivity predictions.) 

Since the column density of hydrogen gas in molecular clouds is well-established from radio observations, 
the uncertainty in estimating the flux of diffuse gamma rays is dominated by the uncertainty in the flux of 
hadronic cosmic rays within the cloud. We can reduce the effect of this uncertainty by considering the ratio 
of gamma ray to cosmic ray fluxes. In the case of a uniform flux of cosmic rays, this ratio will simply be 
proportional to the average column density. For example, using the value of T from Aharonian (1991), the 
expected diffuse emission is given as: 



■f ~2x 10- 5 - 



<N H > 



22- 



(3) 



where < Njj >22=< Nh > /(10 22 cm 2 ) is the average column density of gas in any direction. In the 
plane of the outer galaxy visible by CASA-MIA (40° < I < 200°), < N H >22 varies between 0.8 and 1.2 for 
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|6| < 5° (e.g., Blocmcn ct al. 1984). In this case, the flux ratio, J 7 /J C r, depends only weakly on the exact 
form of the cosmic ray spectrum, and is expected to be constant as a function of energy to within a factor of 
2 over the range from 100 to 400 TeV. A realistic estimate of the emission based upon a three-dimensional 
model of galactic column density yields a similar result ( Berezinsky et al. 1993 ). We note that if the 
ambient flux of cosmic rays within a particular cloud is stronger than the average flux measured at earth 
(e.g., because the cloud is in close proximity to a source of cosmic rays), then we would expect the diffuse 
gamma-ray flux to be stronger than predicted. 



2. Observations and Analysis Method 

When ultrahigh energy cosmic rays and gamma rays strike the upper atmosphere of the earth, they 
produce large cascades of electromagnetic and hadronic particles. For primary particle energies around 100 
TeV these particle cascades, called extensive air showers, reach the ground and can be detected readily 
by surface scintillation counters. At ground level, the showers consist largely of electrons, photons, and 
muons. A vertically incident 100 TeV proton primary produces a shower at ground level of roughly 
30,000 electrons and 2000 muons. The shower consists of a thin "pancake" of relativistic particles, whose 
orientation preserves the direction of the incident primary particle. A sparse array of detectors sampling a 
small fraction of the shower particles can thereby determine the total number of particles in the shower and 
the arrival direction of the primary particle with reasonable precision. The combination of the Chicago Air 
Shower Array (CASA) and the Michigan Muon Array (MIA) comprise the largest and most sensitive such 
array built to date. 

The CASA experiment consists of 1089 scintillation detectors laid out on a square grid of 15m spacing. 
The experiment is located at Dugway, Utah (40.2° N, 112.8° W) at an atmospheric depth of 870 g/cm 2 . 
The total collection area of CASA is roughly 230,000 m 2 , of which 0.7% is actually covered by scintillator. 
The detectors have local analog and digital electronics to record the number of particles and the shower 
arrival time at each scintillator. The number of particles detected can be used to estimate the energy of 
the primary particle (cosmic ray or gamma ray). By comparing the observed shower rate to the known 
cosmic ray flux, we find that the median energy of gamma rays creating showers detectable by CASA is 
approximately 120 TeV QBorionc ct al. 1997b| ). 



At the Dugway site, beneath CASA, we have constructed a very large array of muon counters - 
the Michigan Muon Array (MIA). This array consists of 1024 scintillation counters buried three meters 
below ground level. The counters are grouped into sixteen "patches" and have a total scintillator area of 
2500 m 2 . The amount of earth overburden above MIA results in a muon energy threshold of ~ 0.8 GeV; 
electromagnetic punch-through is negligible for all but the very largest showers. A complete description of 



the C ASA-MI A instrument can be found elsewhere (Borione et al. 1994) 



In showers initiated by cosmic ray nuclei, muons are produced from charged pion and kaon decays 
in the central hadronic core. Showers initiated by gamma rays do not produce muons in significant 
quantities because the cross section for 7-air hadroproduction is much smaller than the cross section for 
electron-positron pair production (7-air — * e + e~). In other words, in contrast to cosmic ray nuclei, gamma 
rays usually interact electromagnetically, and produce air showers with fewer charged mesons, and therefore 
fewer muons at ground level. On average, CASA-MIA detects approximately nine muons for each cosmic 
ray air shower; in contrast, we expect to detect only 0.25 muons for each gamma-ray shower, on average. 



To search for diffuse concentrations of cosmic gamma rays, we employ the technique developed by 
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Matthews et al. (1991). The number of muons resulting from a gamma-ray initiated shower should be, on 
average, much less (by a factor of 30 to 40) than the number of muons from a hadron-initiated shower of 
comparable energy. We therefore search for gamma-ray emission by looking for a localized excess of cosmic 
rays which are muon-poor with respect to the number of muons expected from typical hadron-initiated 



showers (Govault et al. 1991) 



We parameterize the muon content of each shower by the quantity: 

(?V-exp) 

where ^ M _ bs an d ?V_exp represent the numbers of muons actually observed and the number of muons 
expected for hadronic showers. We determine the expected number of muons on a shower-by-shower basis 



by parameterizing the lateral density of muons in terms of the function p M (-R, N e , 9 Z ) described by Greisen 



(1960), where R is the distance to the shower core, N e is the electron shower size, and 9 Z is the zenith angle 



of the incident shower. 

Having parameterized the muon content in this manner, we consider the expected distribution of r M for 
both hadronic and gamma-ray primaries. For hadrons, we expect that the distribution of is characterized 
by the actual array data, since the vast majority of air showers (at least 99.9%) are known to be hadronic. 
For gamma rays we must invoke Monte Carlo simulations to determine the expected distribution of r^. The 
distribution for gamma rays will be shifted significantly towards lower values of r M and will include a much 
larger fraction of events with zero muons detected. The ry, distributions for both hadronic and simulated 
gamma-ray air showers are shown in Figure [lj 

We define as muon-poor those events which have a value of r M less than a value chosen to optimize 
sensitivity to gamma rays. This sensitivity scales as h^/\/h CR , where /i 7 and /i CR represent the fraction 
of gamma rays and cosmic rays retained in the muon-poor sample, respectively. The optimal value of ry, 



depends weakly on the source declination and the energy of the selected events (Borione et al. 1997a) 



To determine if there is evidence for gamma-ray emission from a localized candidate source region of 
the sky, we examine the distribution of r M for showers within this region, and we search for an excess of 
muon-poor events. The excess is determined relative to a comparison region which presumably contains a 
background of uniformly hadronic showers. To compare the r M distributions from source and background 
regions, we normalize so as to match the total number of events with r M > 0.0, assuming that all such 
events are purely hadronic. We note that this technique has no sensitivity whatsoever to the presence of 
any isotropic flux of gamma rays which may be present in our data. 

To the extent that our parameterization of the muon lateral distribution applies generally to all 
showers, we expect that the distribution of r M should remain identical for all hadronic showers at a fixed 
overburden throughout the atmosphere. However, since the overburden changes rapidly with zenith angle 
to the shower, we expect that the shape of the r M distribution will depend upon the zenith angle. If 
uncompensated, these systematic changes in the shape of the r M distribution can have a significant impact 
upon our estimate of muon-poor showers. 

We have developed a technique for characterizing and removing the systematics in the r M distribution. 
We divide the entire data sample into a series of coarse bins in time (about one bin per day) and horizon 
coordinates (8-10 bins each in elevation and azimuth). For each source event, we generate comparison 
background events by sampling a value of r M from the distribution corresponding to the appropriate horizon 
coordinate bin. Thus, systematic changes to the shape of the distribution that are local to the events 
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Fig. 1. — The 77, distribution of simulated gamma ray (dashed) vs. detected hadronic (solid) cosmic rays. 
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in each sky bin will appear in both source and background distributions. This technique removes most 
systematic effects, provided that the bin size is chosen so that the r M distribution obtained from the data 



represents a good approximation of the ry, distribution for the entire bin (Covault et al. 1994) 



Figure || graphically demonstrates the application of this technique to real data. The distribution 
for candidate diffuse source events within the galactic plane (± 5°) is plotted against the distribution 
for the off-source background region. Analysis of residuals shows no major differences between the source 
and background distributions, and gives us confidence that the systematics have largely been removed. 
We attribute the excess in the chi-squared (xt =1.25, v = 51 d.o.f) to either an unremoved systematic 
effect, or the marginal presence of a possible signal. If there were a detectable flux of gamma rays from this 
region of the galactic plane, Figure ^| would have shown an excess of signal events over background on the 
left-hand side of the plot. Indeed, there is an excess (1.63 standard deviations in the plot shown) which, 
while consistent with the expected flux, is not significant enough to claim as a detection. 



3. Results 

We apply this analysis technique to search for evidence of gamma rays from the galactic plane region in 
the energy range from 140 to 1300 TeV. No significant evidence for diffuse emission is found. Upper limits 
to diffuse flux from the plane of the galaxy are calculated taking into account the slight loss in sensitivity 
(< 10%) due to the finite angular resolution of CAS A. We present 90% confidence level upper limits on the 
flux for emission for a range of emission disk thicknesses {\bn\ < 2°, 5°, and 10°) for comparison to previous 



predictions and experiments. (For example, Aglietta et al. 1992 have claimed that the band \bn\ < 2 



best matches the expected UHE emission, based on COS-B results, while we believe that the EGRET 
results demonstrate that \bji\ < 5° is the most appropriate search region). The CASA-MIA flux limits are 
plotted in Figure |^. Also shown are upper limits from previous experiments and a theoretical prediction 
for the expected flux. We note that cosmic ray rejection improves rapidly with energy, thus compensating 
for reduced statistics, so that the sensitivity of CASA-MIA to diffuse emission is relatively constant in the 
region from 200 to 1000 TeV. 



4. Discussion 

The observed weak excess of muon-poor events in several of the selected spatial and energy ranges 
could be taken as evidence for the presence of diffuse emission from the galactic plane. Indeed, if the excess 
is interpreted as a positive gamma-ray signal, the inferred flux is consistent with current model predictions. 
However, in view of the uncertainty in assessing the impact of any remaining systematics, we interpret 
the observed excesses as not statistically significant enough to warrant a claim for detection, and we use 
our results to place upper limits on the predicted emission. We note that future prospects for increased 
sensitivity to verify or refute the presence of a possible signal at this level are not good, since this would 
require a much larger experiment and/or a much longer running time, neither of which is foreseen. 

The CASA-MIA flux limits place strong constraints on models for diffuse emission. These limits are 
below the predictions of some earlier models (e.g., Bcrczinsky and Kudryavtsev 1990| ) and are approximately 



a factor of two above the minimum emission predicted by more recent models (e.g., Aharonian 1991 



Berezinsky et al. 1993). These predictions are based upon diffuse emission due to the process of normal 



cosmic ray nucleon-nucleon interactions on passive gas targets only. Therefore we conclude that no other 
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Fig. 2. The upper plot shows the muon content (r^ distribution) for CASA-MIA air shower events 
arriving from the direction of the galactic plane (\bu\ < ±5°). Points with error bars represent the on- 
source distribution, while histogram lines represent the off-source distribution. The lower plot indicates the 
on-source minus off-source residuals (sigmas per bin). Evidence for the presence of ultra- high energy gamma 
rays would be indicated by an excess of muon-poor events. A statistically insignificant excess (1.63er) is seen. 



- 8 - 



i-2 



Limits (90% C.L.) to Diffuse y— Rays from Plane 



t 1 1 — r 



,-3 



-5 



BASJE 1991 



GREX 1995 



EAS-TOP 1992 



t UM Array, 1991 

. CASA= MIA 1996 



I 



-> 



Predicted .diffuse flux 
for galactic plane ± 5 degrees 
as seen from Northern Hemisphere 



J I I I I L 



2 



Energy (TeV) 



Fig. 3. — CASA-MIA sensitivity to diffuse gamma ray emission from the central plane of the galaxy 
(\bu\ < ±5°, 50° < £u < 200°). Sensitivities are given in terms of the fraction of gamma rays relative 
to the detected all-particle flux of cosmic rays at the earth. Also shown are limits from previous experiments 
(BASJE-[Kakimoto et al. 199l|, EAS-TOP-|Aglietta et al. 1992], UM-|Matthews et al. 1991].) Predicted flux 



from Aharonian 1991 
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Table 1. Limits to Diffuse Emission 



Region (50° < £ u < 200°) Median Energy (TeV) Significance (cr) J 7 /J C R 90% C.L. (1CT 5 ) 



-2° < bn < 2° 


140 


+ 1.78 


7.2 




180 


+ 1.81 


3.8 




310 


+2.56 


5.2 




650 


+ 1.12 


3.2 




1300 


+0.07 


4.6 


-5° < bn < 5° 


140 


+ 1.63 


3.4 




180 


+0.08 


2.6 




310 


+0.86 


2.4 




650 


+ 1.60 


2.6 




1300 


+0.06 


3.5 


-10° < b u < 10° 


140 


+2.39 


2.8 




180 


+ 1.79 


2.2 




310 


+0.87 


2.3 




650 


+0.91 


1.8 




1300 


-0.56 


2.3 



Note. - Tabulated upper limits to diffuse gamma-ray emission from the plane of the galaxy. 
Although positive excesses are seen, we do not view these as statistically significant enough to claim 
detections. Flux limits are tabulated for bands along the galactic plane from \b u \ < 2° to |6 n | < 10°. 
Median energy is quoted for integral flux limits. Selected spatial regions and energy bands are not 
statistically independent. 
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emission process is likely to dominate at ultrahigh energies. Some models for cosmic ray acceleration 
suggest that gamma rays from clouds in proximity to SNR may be enhanced and/or spectrally hardened 
(e.g., Drury, Aharonian, and Volk 1994). Others have suggested that inverse Compton scattering may 



be a significant component of galactic diffuse emission at ultrahigh energies (e.g., Aharonian and Atoyan 
1995| ). A report by Bloemen 1991 based upon an analysis of the COS-B data suggests evidence for a 



spectral hardening at energies above 1 GeV (A 7 = 0.4) for the outer galaxy; however this result has not 
been confirmed by EGRET ( [Hunter et al. 1997| ). The CASA-MIA upper limits rule out any such spectral 
hardening above 200 TeV in the outer galaxy with A 7 > 0.1. 
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